Is There a Higher-Order Mode Coupling Transition in Polymer Blends? 
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We present simulations on a binary blend of bead-spring polymer chains. The introduction of 
monomer size disparity yields very different relaxation times for each component of the blend. 
Competition between two different arrest mechanisms, namely bulk-like dynamics and confinement, 
leads to an anomalous relaxation scenario for the fast component, characterized by sublinear time 
dependence for mean squared displacements, or logarithmic decay and convex-to-concave crossover 
for density-density correlators. These anomalous dynamic features, which are observed over time 
intervals extending up to four decades, strongly resemble predictions of Mode Coupling Theory for 
nearby higher-order transitions. Chain connectivity extends anomalous relaxation over a wide range 
of blend compositions. 

PACS numbers: 64.70.Pf, 83.80.Tc, 83.10.Rs 



I. INTRODUCTION 

The widespread potential applications of polymer 
blends have spurred a great deal of recent research 
on the dynamics of these systems showing tunable 
properties. Dynamics in polymer blends are directly 
related with one fundamental question of polymer- 
based soft matter systems: the way local friction 
arises in complex environments that are chemically 
heterogeneous. Thermodynamically miscible polymer 
blends exhibit complex dynamics. They show 'dy- 
namic heterogeneity' in the meaning that two different 
segmental dynamics (a-relaxation) can be observed 
in the blend. Moreover, a symmetric broadening of 
any relaxation function of each component is usually 
observed as the glass transition of the blend is ap- 
proached. These features are usually understood in 
terms of self-concentration effects induced by chain 
connectivity ll| and thermally driven concentration 
fluctuations [2|. 

For most of the investigated systems, the two com- 
ponents in the blend display qualitatively similar dy- 
namic features. However, recent experimental results 
by nuclear magnetic resonance Q and neutron scatter- 
ing suggest that a rather different picture emerges 
when the two blend components exhibit very different 
segmental mobilities in the homopolymer state and 
the concentration of the fast component in the blend 
is rather dilute. Then, the segmental dynamics of the 
fast component strongly deviate from the expected 
behavior and seem to be qualitatively different from 
that of the slow component. Moreover, fully atomistic 
molecular dynamics simulations (MDS) suggest the 
possibility of an unusual logarithmic time decay for 
the intermediate scattering function of the fast com- 
ponent 0. 

Simulations are a valuable tool for investigating the 
microscopic origin of these unusual features. With 
present computational resources, a systematic study 
— i.e, by exploring a wide range of control parame- 
ters as temperature or blend composition — can only 
be performed by using very simple models, providing 
they display the essential dynamic features observed 



in real systems. With these ideas in mind, in this ar- 
ticle we present an extensive MDS investigation on a 
blend of bead-spring polymer chains with components 
of very different mobilities. The fast component ex- 
hibits an anomalous relaxation scenario rather differ- 
ent from that characteristic of usual liquid-glass tran- 
sitions. Namely, mean squared displacements display 
sublinear behavior, and density-density correlators ex- 
hibit a convex-to-concave crossover, which is observed 
both by varying wavevectors and control parameters. 
By properly tuning the latters, a logarithmic decay 
can be obtained. 

These anomalous features are observed at interme- 
diate time intervals extending up to four decades, and 
display striking analogies with Mode Coupling Theory 
(MCT) predictions for states close to a higher-order 
MCT transition, which emerges as the result from the 
competition between different arrest mechanisms - 
we suggest bulk-like caging and confinement for the 
present case. 

The article is organized as follows. In Section II 
we introduce the model and give computational de- 
tails. In Section III we present simulation results and 
evidence that the fast component displays unusual re- 
laxation features. In Section IV the framework of the 
Mode Coupling Theory is used in an operational way 
to describe simulation results. Conclusions are given 
in Section V. 

II. MODEL AND SIMULATION DETAILS 

We have carried out simulations on a model very 
similar to one originally introduced by Grest and Kre- 
mer |5j for homopolymers. Here we investigate a bi- 
nary mixture of bead-spring polymer chains. Each 
chain consists of N = 10 identical monomers of the 
species A or B. Monomers of both species have equal 
mass m = 1. The interaction between two given 
monomers, placed on different chains or in a same one, 
is given by a soft-sphere potential plus a quadratic 
term: 



V„r = Ae\(a n alr) 
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where e = 1 and a, (3 s {A, B}. The interaction is 
zero beyond a cutoff distance ca a {3, with c = 1.15. 
The vaiues C — 7c -12 and C 2 = 6c~ 14 guarantee 
continuity of potential and forces at the cutoff dis- 
tance. The radii of the soft potential for the different 
types of interaction are <tbb = 1, oaa = <5<7bb, and 
cab = (caa + 0bb)/2. The potential is purely 
repulsive. It does not show local minima within the 
interaction range r < co a @. Hence, slow dynamics in 
the present model arises as a consequence of steric ef- 
fects, as in the computationally more involved case of 
chains of hard-sphere monomers. 

Adjacent monomers of a same chain also interact 
through a^finite extensible nonlinear elastic potential 
(FENE) " 
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FENE 



(r) = -A;J^eln[l-( J R a ac( )-V], (2) 



where k = 15 and Rq = 1.5. The superposition of the 
potentials QJ and provides an effective bond po- 
tential for adjacent monomers with a sharp minimum 
at r = 0.985<7 QQ , which makes bond crossing impos- 
sible for the investigated values of the monomer size 
disparity 5. 

The blend composition is defined as xb = 
Nb/(Na+Nb), with Aa and Ab denoting respectively 
the number of A- and B-chains. Data here presented 
for different compositions correspond to a single value 

5 = 1.6. Simulations were also carried out for values 

1 < 8 < 1.75 at a single composition ib = 0.5. For 

6 < 1.2 we did not observe signatures of the anoma- 
lous relaxation features exposed below. 

The packing fraction, (f>, is defined as <j) = 
(■K/6)L~ 3 [N A a 3 AA + N B a 3 BB ], with L the side of the 
simulation box. In the following, temperature, T, dis- 
tance and time, t, will be given respectively in units of 
e/fcg, obb and (TB^m/e) 1 / 2 . Simulations have been 
carried out at a constant packing fraction <f> = 0.53. 
This value is comparable to those used in simulations 
of slow relaxation in simple liquids Q • 

We investigate the T-dependence of the dynamics 
for blend compositions in the range 0.1 < xb < 1. 
The total number of chains at each composition varies 
between 250 and 600, with a minimum of 60 for the 
minoritary component. The system is prepared by 
placing the chains randomly in the simulation box, 
with a constraint that avoids core overlapping. Pe- 
riodic boundary conditions are implemented. Equa- 
tions of motion are integrated by using the veloc- 
ity Verlet scheme with a time step ranging from 

2 x 10~ 4 to 5 x 10~ 3 , for respectively the highest and 
the lowest temperature. A link-cell method [S| is used 
for saving computational time in the determination of 
the particles within the interaction range of a given 
one. 

At each state point, the system is thermalized at 
the requested temperature by periodic velocity rescal- 
ing. After reaching equilibrium, energy, pressure, 
chain radii of gyration and end-to-end distances show 
no drift. Likewise, dynamic correlators and mean 
squared displacements show no aging, i.e., no time 
shift when being evaluated for progressively longer 
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FIG. 1: Symbols: mean squared displacement of B- 
particles for the homopolymer (a) and for a blend compo- 
sition xb = 0.3 (b). Thick straight lines correspond to the 
linear diffusive regime ((Are) 2 ) oc t. Thin straight lines 
in panel (a) correspond to sublinear behavior ((Are) 2 ) oc 
t ' 65 . Thin straight lines in panel (b) correspond to 
((Are) 2 ) cx t a with a T-dependent exponent a. From top 
to down, a = 0.635, 0.58, 0.47, 0.44 and 0.415. Results 
for A-particles (curves) at xb = 0.3 are also shown for 
comparison. The exponent for the observed intermediate 
sublinear regime is 0.66. 



time origins. Once the system is equilibrated, a mi- 
crocanonical run is performed for production of con- 
figurations, from which dynamic correlators and mean 
squared displacements are computed. For each state 
point, the latter quantities are averaged over typically 
20-40 independent samples. 

III. SIMULATION RESULTS 

a) Mean squared displacements 

Fig. shows the T-dependence of the mean 

squared displacement of the B-particles, ((ArB) 2 }, for 
the homopolymer case (xb = 1)- We observe fea- 
tures analogous to similar bead-spring models previ- 
ously investigated . A plateau is developed after the 
initial ballistic regime (t < 0.2), extending over longer 
times with decreasing temperature. This plateau cor- 
responds to the well-known caging regime — i.e., the 
temporary trapping of each particle induced by their 
neighbors — observed when approaching a liquid-glass 
transition [9j, llCl . At longer times a crossover to a 
Rouse-like 0, UM sublinear regime ((Ar B ) 2 ) oc i 65 
occurs, before the final crossover to linear diffusion at 
very long times, which is reached only for the highest 
investigated temperatures. The T-independent expo- 
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nent 0.65 is compatible with previous simulations of 
similar models of bead-spring homopolymers [|J . 

Fig. ^ shows results for A- and B-particles for a 
blend composition xb = 0.3. Monomer size disparity 
introduces a clear separation in the time scale of both 
components. Hence, on average, the big A-particles 
and the small B-ones are respectively slow and fast 
particles. A-particles exhibit features analogous to 
those observed in the homopolymer case. It is note- 
worthy that this observation is fulfilled in the whole 
composition range, suggesting that bulk-like dynamics 
is the only relevant arrest mechanism for the slow com- 
ponent. However, in the range xb ^ 0.8, B-particles 
exhibit a rather different behavior, as illustrated in 
Fig. ^ f° r a composition xb — 0.3. The arrest mecha- 
nism is characterized by a larger localization length as 
compared to the homopolymer or to A-particles. This 
is evidenced by the larger value of the mean squared 
displacement at the inflection point after the end of 
the ballistic regime. A direct crossover from the bal- 
listic regime to sublinear behavior ((Are) 2 ) cx t a oc- 
curs, but differently from the homopolymer or the A- 
particles, the exponent a is not constant, but shows a 
clear decrease with decreasing T. At sufficiently high 
temperature, again a ~ 0.6, — i.e., above a certain T 
there are no essential differences between the trans- 
port mechanism for both species. However, for the 
lowest investigated temperatures we find a ~ 0.4. 

b) Density-density correlators 

We compute partial density-density correlators 
F a p(q,t) = (p a (q,t)pp(-q,O)}/(pa,(q,O)p0(-q,O)), 
with p a (q, t) — Sj exp[iq • r QJ -(i)], the sum extending 
over all the particles of the species a e {A,B}. Re- 
markable differences between the homopolymer and 
the blend are observed for density-density correlators 
for B-B pairs, Fbb(<Z, i)- Figs. and[2t> show sim- 
ulation results as a function of T respectively for the 
homopolymer and for a blend composition xb = 0.3. 
The wavevector q corresponds to the maximum of the 
total static structure factor. Figs. OK and Ob show 
the corresponding g-dependence at two selected low 
values of T. The homopolymer exhibits the standard 
behavior in the proximity of liquid-glass transitions 
H H El EH E3- After the initial transient regime, 
Fbb(<Z, t) shows a first decay to a plateau. With de- 
creasing temperature, the plateau extends over longer 
time intervals. At long times, a second decay occurs 
from the plateau to zero. This second decay corre- 
sponds to the a-process of the glass transition and 
is well described by a stretched exponential function. 
Analogous results are obtained for Faa (q, t) in the 
whole composition range, as expected from the results 
above presented for the mean squared displacement of 
A-particles. 

Rather different features are observed for Fbb (<Z, t) 
for compositions xb < 0.8. As can be seen in Fig. 
Eb, Fbb does not exhibit a defined plateau at 
any value of q or T . Moreover, an unusual convex- 
to-concave crossover, both in q and T, is observed 
for the shape of the correlators, and in particular 
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FIG. 2: Symbols in panels (a) and (b): T-dependence 
of FBs{q,t), respectively for the homopolymer and for a 
blend composition xb = 0.3, as computed from the simula- 
tions. The straight line in panel (b) indicates logarithmic 
decay over four time decades. For xb = 0.3, the correla- 
tor FA\{q,t) is also shown (curves) for comparison. The 
wavevectors are q — 7.1 for the homopolymer and q = 4.5 
for xb = 0.3, both values corresponding to the maxima of 
the respective total static structure factors. 

there are state point-dependent values of q for which 
logarithmic relaxation occurs over time intervals ex- 
tending up to four decades. Recent atomistic simu- 
lations Q, at a single temperature and composition, 
of (fast) poly(ethylene oxide), PEO, in a (slow) ma- 
trix of poly(methyl methacrylate) report an impres- 
sively similar crossover at the g-dependence for self- 
correlations of PEO hydrogens (see Fig. 8 in Ref . 4) , 
supporting the simple bead-spring model for polymer 
blends here investigated. 

IV. DISCUSSION 

The unusual features observed for mean squared 
displacements and dynamic correlators of the fast 
component strongly resemble those reported by sim- 
ulations [H El EH and experiments El El El 
for hard-sphere colloids with short-ranged attrac- 
tive interactions. For these systems, such features 
have been recently rationalized in terms of the MCT 
El El HI El H|. Motivated by this fact, we 
discuss the present results by using MCT in an oper- 
ational way. In its ideal version — i.e., not inclu ding 
activated hopping events — MCT predicts H 0, Ea | 
a sharp transition from an ergodic liquid to a non- 
ergodic arrested state ('glass') at a given value x c of 
the relevant control parameter x (in practice (f> or T) . 
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FIG. 3: Symbols in panels (a) and (b): q-dependence 
of FBs{q, t), respectively for the homopolymer and for a 
blend composition ib = 0.3, as computed from the simu- 
lations at two selected low temperatures. Curves in panels 
(a) and (b) are respectively fits to Eqs. ||3J and (J1J. 

When crossing the transition point from the ergodic to 
the arrested state, the long-time limit of the density- 
density correlator for wavevector q, F(q, t), jumps dis- 
continuously from zero to a finite value, denoted as the 
critical non-ergodicity parameter, f°. In the MCT for- 
malism, the standard liquid-glass transition is of the 
fold type (also denoted as A 2 ) @, ES IH El- The 
initial part of the a-process (von Schweidler regime) 
is approximated by a power law expansion [9J, llfj . |l2( : 

F{q,t)^r q -h q {tM h + hf){t/rf\ (3) 

(2) 

where the prefactors h q and h q only depend on the 
wavevector q. The characteristic time r only depends 
on the separation parameter \x — x c \ and diverges at 
the transition point. The exponent b is related to the 
so-called exponent parameter 1/2 < A < 1, via the 
relation A = r 2 (l_ +- b)/T(l + 2b), with T the Gamma 
function HUSO!- 

As previously found for similar bead-spring models 
[H EH- Eq. © provides a good description for the 
decay from the plateau in the homopolymer case (see 
Fig. [3^). A value b — 0.477 is obtained. The corre- 
sponding exponent parameter is A = 0.799. The value 
of the constant exponent a = 0.65 for the observed 
Rouse-like regime in the mean squared displacement 
is consistent with MCT theoretical calculations for 
chains of N — 10 hard-sphere monomers (a = 0.6) 
|25) . The obtained q- modulation for is also consis- 
tent with MCT expectations and previous inves- 
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FIG. 4: For the homopolymer (xb = 1), temperature de- 
pendence of the a-relaxation time of FsB(Q,t) at q corre- 
sponding to the maximum of the static structure factor. 
Symbols are simulation data. The line is a fit to a power 
law to. 2 oc (T — T c ) -7 (see text). 

tigations of bead-spring models E3 ■ 

Fig. 0] shows, for the homopolymer, the In- 
dependence of the relaxation time of Feb (q,t), for q at 
the maximum of the static structure factor. We define 
the former as the time to. 2 where Fbb (o, t) decays to 
0.2, a value well below the plateau. This corresponds 
to times probing the a-relaxation. According to MCT 
for A 2 -transitions 0, 0, 0], such times diverge at 
the transition point as a power law oc \x — x c |~ 7 (in 
the present case x = T). The exponent 7 is related 
to A through the expressions 7 = 1 /2a + 1/26 and 
A = T 2 (l - a)/r(l - 2a). From the value A = 0.799 
given above we obtain 7 = 2.84. Consistently, a good 
description of the data is obtained by forcing this value 
of 7 and leaving T c as a free parameter (see Fig. 0J). 
We obtain T c = 0.517. 

The features observed for B-particles in the compo- 
sition range xb J$ 0.8 cannot be rationalized within 
the picture of a MCT A2-transition. However, they 
exhibit a striking resemblance with MCT predictions 
in the proximity of higher-order transitions A„ + i, 
characterized by A = 1 , which can emerge as the result 
from the interplay between n > 2 control parameters 
x±, xi ... x n . Higher-order MCT transitions were ini- 
tially derived for schematic models I27l l2 8l. an d later 
for short-ranged attractive colloids |l9j, |2JJ, |2lJ, |22j, |2jj 
as a first realization in real systems. Close to a higher- 
order transition, or more generally to a fold transi- 
tion with A < 1, an anomalous relaxation scenario 
emerges. The mean squared displacement exhibits an 
intermediate sublinear regime ji^, with a decreas- 
ing exponent as the transition is approached. As re- 
ported above, this behavior is observed for .B-particles 
for blend compositions xb < 0.8. In the higher- 
order MCT scenario, F(q, t) does not exhibit a defined 
plateau. Instead, in an intermediate time interval of 
several decades^it is approximated by a logarithmic 
expansion [H |H El IH3 : 

F(q, t) » f c q - H q ln(i/r) + Hf \n 2 (t/r), (4) 

(2) 

where the prefactors H q and Hq depend on q and 
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FIG. 5: Panels (a), (b), and (c) are respectively the values 
of fq, H q , and H q 2 \ for two different state points SI and 
S2, as obtained from fits of F BB (q, t) to Eq. Q). The state 
points are Si: (x B = 0.3, T = 0.4) and S2: (x B = 0.3; 
T — 0.5). The characteristic times r for SI and S2 are 
respectively 1100 and 67. The fq for the homopolymer 
(xb = 1), as obtained from fitting to Eq. is also shown 
for comparison in panel (a). The dashed line in panel (c) 
indicates zero value for H q 2 ^ . 



on the distance of the state point {x 71 } to the transi- 
tion point {x™}. As shown in Fig. 03, Eq. 10} pro- 
vides an excellent description of the simulation data 
for -Fbb(<7, t) in the blend. It is worthy of remark that 
the observed convex-to-concave crossover is present in 
the higher-order MCT scenario dHUHl]. It is in- 
deed one of its main signatures and differentiates it 
from other theoretical frameworks ■ The pure log- 
arithmic decay observed at some T- and q-values is 
also present in the higher-order MCT scenario, which 
predicts [51 EtI Ejf lines in the control parameter 

space with H q ' =0. 

Figs. |S] and [B] show the values of f q , H q , and Hq 
obtained from the fits to Eq. J3J of Fbb(<7, i) for two 
sets of state points [S1,S2] and [S3,S4]. The coordi- 
nates of such state points are S1:(ccb = 0.3, T — 0.4), 
S2:{x B = 0.3, T = 0.5), S3:(x B = 0.6, T = 0.4), 
and S4:(a;B = 0.6, T = 0.6). The fact that a com- 
mon value of /| is found for the states of each set, 
together with the scaling behavior of H q (see below) 
would be consistent with the existence, for each set, 
of a nearby MCT higher-order transition (A = 1) or 
a fold transition with A < 1. Interestingly, f q shows 
no strong modulation but a nearly monotonous de- 
cay (see Figs. andEt): resembling the qualitative 
behavior observed at A3- and A4-transitions in short- 



FIG. 6: As Fig. Q>] for the state points S3: (x B = 0.6, 
T = 0.4) and S4: (x B = 0.6; T = 0.6). The characteristic 
times t for S3 and S4 are respectively 200 and 7.4. 

ranged attractive colloids El El El HI El • 

According to MCT ElEJ the coefficient H q in Eq. 
J3J factorizes as H q — C({x n })H q , where H q only 
depends on q, and the q-independent term C({x"}) 
depends on the state point. Hence, the values of H q 
obtained for different state points close to the transi- 
tion point must be proportional. This scaling behav- 
ior is indeed displayed by simulation data, as demon- 
strated in Figs. 03 andlHb. As for f q , the qualita- 
tive q-dependence of H q resembles that observed at 
higher-order MCT transitions in short-ranged attrac- 
tive colloids. 

Also resembling MCT predictions for the higher- 
order scenario 0, l^l , the coefficient H q 2 ^ is smaller 

than H q , and does not obey scaling. Indeed any test 

(2) 

of scaling behavior over a wide q-range for Hf' is un- 
satisfactory (see two examples in Figs. |SJ; andEt). 

(2) 

Apparently, the value of q for which H q ' = is de- 
pendent on the state point, in agreement with MCT 
expectations. 

The striking dynamic analogies of the bead-spring 
polymer blend here investigated with the higher-order 
MCT scenario, which has been derived for short- 
ranged attractive colloids, provide connections with 
different research fields. In these latter systems, in a 
narrow range of cj> and T, anomalous relaxation origi- 
nates from the competition between two arrest mecha- 
nisms of very different localization lengths: reversible 
bond formation induced by the short-ranged attrac- 
tive potential, and hard-sphere repulsion [3 E El 
Elj| . The dynamic analogies between so different sys- 
tems suggest that this anomalous scenario might be 
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a general feature of systems with several competing 
mechanisms for dynamic arrest, and should motivate 
theoretical approaches within the MCT framework. In 
the case here investigated, we suggest bulk-like caging 
and confinement as the two competing mechanisms for 
arrest of the fast B-component. While the former is 
induced by the neighboring B-particles, the latter is 
induced by the matrix formed by the chains of the 
slow A-component. 

It is worthy of remark that, due to chain connectiv- 
ity, any given particle has always neighboring particles 
of its same species. Hence, competition between bulk- 
like caging and confinement for the dynamics of the 
fast component occurs even at high dilution of the B- 
particles. Indeed, all the reported anomalous features 
are observed in a composition range from xb ^ 0.8 to, 
at least, the smallest investigated value xb =0.1. 

Results reported here must not be understood as a 
proof of the existence of higher-order MCT transitions 
(or fold transitions with exponent parameter A < 1) in 
the system here investigated. A proper answer to this 
question could only be provided by solving the corre- 
sponding MCT equations. Still, it is worth mention- 
ing that recent theoretical calculations by Krakoviack 
32] have reported a higher-order MCT transition for 
a binary mixture of mobile and static hard spheres, 
which can be seen as a simplified model for liquids in 
confining media with interconnected voids. This sys- 
tem, where competition between bulk-like dynamics 
and confinement occurs, is a priori closer to polymer 
blends than short-ranged attractive colloids. This re- 
sult supports the hypothesis of an underlying higher- 
order MCT scenario in the model here investigated. 

V. CONCLUSIONS 

We have carried out simulations on a bead-spring 
model for polymer blends. This model displays non- 
conventional glassy dynamics. Mean squared displace- 
ments and density-density correlators for the fast com- 
ponent exhibit many of the precursor effects char- 
acterizing nearby MCT higher-order transitions, or 
fold transitions with exponent parameter very close 
to unity. 

These anomalous features are analogous to those 
observed for short-ranged attractive colloids, for 
which a higher-order MCT scenario has been de- 
rived, suggesting that such features might be gener- 
ally present in systems with several competing mecha- 
nisms for dynamic arrest. We suggest bulk-like caging 
and confinement for the case of polymer blends. Chain 
connectivity extends competition between both mech- 
anisms over a wide range of blend compositions. The 
strong resemblance of dynamic correlators with those 
recently reported in fully atomistic simulations of a 
real polymer blend strongly supports the dynamic pic- 
ture here presented. These results open the possibility 
of new theoretical approaches for investigating poly- 
mer blend dynamics. 
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